Universal doping dependence of the ground state staggered magnetization in cuprates 
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Using muon spin rotation we determine the zero temperature staggered antiferromag- 
netic order parameter Mo versus hole doping measured from optimum Ap m , in the 
(Ca a; Lai_ a ;)(Bai.7 5 _ a; Lao.25+a:)Cu30j / system. In this system the maximum T c and the superex- 
change J vary by 30% between families (a;). Mo(a;,Ap m ) is found to be ^-independent. Using 
neutron diffraction we also determine the lattice parameters variations for all x and doping. The 
oxygen buckling angle is found to change with x, implying a change in the holes kinetic energy. We 
discuss the surprising insensitivity of Mo(x, Ap m ) to the kinetic energy variations in the framework 
of the t-J model. 



It is widely agreed that the antiferromagnetic (AFM) 
phase of the cuprates should be addressed as a doped 
Mott insulator, where holes are moving on a 2D AFM 
background This scenario is described by the t-J 
model Hamiltonian where t and t' are the near and next- 
near neighbor hoppings, respectively, and J is the Heisen- 
berg superexchange. Above some critical doping the zero 
temperature staggered AFM order parameter Mo is de- 
stroyed, and the cuprates enter a glassy, phase separated, 
state. Since the glassy state precedes superconductivity, 
understanding this transition is crucial to understand- 
ing the cuprates. Particularly interesting is the doping 
dependence of Mq and its variations with the different 
energy scales. These variations were calculated theoreti- 
cally [2j but not measured in a controlled manner. Such 
measurements could shed light on the effective Hamilto- 
nian governing the holes at T — » in the underdoped 
region. While J can be measured relatively simply with 
neutron or Raman scattering on a single crystal, it is very 
difficult to determine t experimentally. 

In this work we determine t/J from their lat- 
tice parameter dependence, including the buckling 
angle, for different cuprate families (x) of the 
(Ca a! Lai_ a! )(Bai.75_ a; Lao.25+a)Cu30 B (CLBLCO) sys- 
tem, where the maximum T c [x/J lax ) and J varies by 
about 30% between families [IJ. This is done by Rietvelt 
refinement of neutron diffraction. We also determine the 
doping dependence of Mq, which is expected to depend 
on t/J. Zero field muon spin rotation (/iSR) is employed 
for this purpose. Our main finding is that the doping 
dependence of Mg is universal despite the fact that t/J 
varies between families. 

We chose to work with the CLBLCO system for sev- 
eral reasons. This is a high temperature superconductor 
(HTSC) with YBa 2 Cu 3 O a (YBCO) structure. The fam- 
ily index x varies in the range 0.1 < x < 0.4. Doping is 
possible all the way to the overdoped regime. All com- 
pounds are tetragonal and there is no chain ordering as in 
YBCO Q. As we show below, there are minimal struc- 
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FIG. 1: (a) The zero temperature muon oscillation angular 
frequency ui as a function of the chemical doping y for all 
four of the (Ca a; Lai_ a; )(Bai.75-^Lao.25+a ; )Cti30 s (CLBLCO) 
families. The antiferromagnetic zero temperature order pa- 
rameter Mo is proportional to ui. (b) The CLBLCO phase 
diagram from [3j including magnetic (close symbols) and su- 
perconducting (open symbols) critical temperatures. The in- 
sets show raw muon polarization data and fits to Eq. [3] for 
three samples from the x — 0.4 family marked in panel b. 



tural differences between the families. In addition, the 
level of disorder as detected by Ca NMR [j| and Cu NQR 
@ is identical for the different families. The phase dia- 
gram is presented in Fig. [ljb) showing the antiferromag- 
netic Neel temperature Tjv, the spin glass temperature 
T g where islands of spins freeze, and the superconduct- 
ing critical temperature T c ; note the axis breaker. In this 
phase diagram Tn Q and T g [?| were measured by /iSR, 
and T c was measured by resistivity jij. The spin glass 
phase penetrates into the superconducting phase. It also 
slightly penetrates into the Neel phase in the sense that 
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a first transition, to long range order, takes place near 
200 K, and a second transition, with additional sponta- 
neous fields, takes place near 10 K. Each transition is a 
continuous function of doping. 

The neutron powder diffraction experiments were per- 
formed at the Special Environment Powder Diffractome- 
ter at Argonne's Intense Pulsed Neutron Source (see Ref 
for more details). Fig. [2] shows a summary of the 
lattice parameters. The empty symbols represent data 
taken from Ref. All the parameters are family- 

dependent, however, not to the same extent. The lattice 
parameters a and c, depicted in Fig.^a) and (b), change 
by up to about 0.5% between the two extreme families 
(x = 0.1 and x = 0.4). The in-plane Cu-O-Cu buckling 
angle is shown in Fig. [2jc). This angle is non-zero since 
the oxygen is slightly out of the Cu plane and closer to the 
Y site of the YBCO structure. The buckling angle shows 
strong variation between the families; there is about a 
30% change from the x = 0.1 family to x = 0.4. This 
change is expected since as x increases, a positive charge 
is moving from the Y to the Ba site of the YBCO struc- 
ture, pulling the oxygen toward the plane and flattening 
the Cu-O-Cu bond. 

We believe that this property is the main cause for 
the different J and therefore different T™ ax between 
the CLBLCO families Nevertheless, we note that 
Pavarini et al. Q showed that in many cuprates families 
rpmax sca j es w ith t' It. t' is controlled by the hybridiza- 
tion of the Cu 4s with the apical oxygen 2p z , hence T™ ax 
scales with the distance R24 between the in-plane copper 
and the apical oxygen. In Fig. [2][d) we show R24 for our 
CLBLCO samples. Our results also support qualitatively 
Pavarini's conclusion. 

From the lattice parameters and the buckling angle it 
is possible to construct the t/J ratio assuming that the 
Hubbard U and the charge transfer energy A are family- 
independent. The basic quantity is the hopping integral 
t p d between a Cu 3d x 2_ y 2 and O 2p [10]. This hopping 
integral is proportional to bond length to the power -3.5 
[HI ]. The hopping from the O 2p to the next Cu 3d x 2_ y 2 
involves again the bond length and cosine of the angle. 
Thus, the Cu to Cu hopping depends on a and 9 as 
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J is proportional to t^ d , hence, 
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Estimates of the tdd and J, normalized to the averaged 
values of the x — 0.1 family, (idd)o 1 and ( J) 1 are pre- 
sented in Fig. [31(a) and (b). Although there is a variation 
in t and J within each family, the variation is much larger 
between the families. J increases with increasing x, in 
qualitative agreement with experimental determination 




FIG. 2: The parameters extracted from a neutron diffraction 
experiment as a function of oxygen doping for the four families 
of CLBLCO. (a) The lattice parameter a. (b) The lattice 
parameter c. (c) 6 - the buckling angle between the copper 
and oxygen in the plane, (d) R24 - The distance between the 
in-plane copper and the apical oxygen. The empty symbols 
are measurements taken from [jj. The lines are guides to the 
eye. 



of J In Fig. [3][c) we show the normalized ratio t/J 
for the four CLBLCO families. There is about 5% dif- 
ference between the two extreme families. But we stress 
that this determination of t/J is only an estimate, used 
in practice to set the oxygen level spacing between sam- 
ples in the /iSR experiment. More accurate calculations 
of t/J are in progress and preliminary data indicate that 
t/J varies by more than 10% between families [T3 ] . 

Next we determine the doping dependence of the order 
parameter using zero field /iSR. The experiments were 
done on the GPS beam line at the Paul Scherrer In- 
stitute, Switzerland. The muon oscillation angular fre- 
quency uj is proportional to the local magnetic field it 
experiences. Therefore, it can be used to determine the 
staggered magnetization M. Typical muon polarization 
curves at T — 5 K are presented in the insets of Fig. [1] 
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FIG. 3: (a) The hopping rate tdd, and (b) the superexchange 
coupling J calculated from Eq. ([!} and ([2}, using the param- 
eters shown in Fig. [2] (c) The ratio tdd/ J as a function of 
doping. The dotted lines are guides to the eye. All data sets 
are normalized to the x = 0.1 familiy. 

for three samples from the x — 0.4 family on the border 
between Neel and glass order. These three samples are 
marked on the phase diagram in Fig. [Tfb). More raw 
data can be found in Ref. [3j. The sample in inset (a) 
is in the spin glass phase; it has no long-range magnetic 
order and hence has no oscillations. The sample in inset 
(c) is in the antiferromagnetic phase, and so it has strong 
oscillations at low temperatures. Finally, the sample in 
inset (b) is an example of an intermediate sample and 
thus has weaker oscillations. 

The best fit of the polarization is achieved with the 
function 

3 

P(t) =Y, A i exp(-A,t) oosfot) (3) 

i=l 

with U13 = 0; the fit is shown in insets (b) and (c) of 
Fig. Q] by the solid line. The reason for multiple fre- 
quencies is that the muons stop at different sites in the 
unit cell. The order parameter extracted from the high 



angular frequency, around a few tens of of radians per 
micro-second, is known to agree with neutron scattering 
experiments 1 1 31 ] . The lower angular frequency is believed 
to emerge from metastable muon sites and is not used for 
further analysis. 

The muon polarization was measured at low temper- 
atures, typically from 5 K to 200 K, and the oscillation 
for T — > was extracted from extrapolation. Fig. [TJa) 
shows a summary of the oscillation angular frequency 
oj(T — > 0) as a function of the chemical doping y for all 
four CLBLCO families. In this plot the AFM critical 
doping, where the oscillations disappear, is different for 
each family. Not surprisingly, this is the same oxygen 
doping where the Neel order is replaced by the spin glass 
phase in the diagram. However, the chemical doping is 
different from the mobile hole doping p m , and a rescaling 
of the doping axis is required. 

The scaling is designed so that the T c domes of all 
families, normalized by T™ ax , will collapse on to a single 
dome. For this purpose the mobile hole parameter mea- 
sured from optimum Ap m is defined by Ap m = K(x) ■ 
(y—ymax), where y ma x is the optimal oxygen doping, and 
K(x) is a family-dependent scaling parameter. K should 
be thought of as doping efficiency parameter connecting 
oxygen level to mobile holes in the CuC>2 planes. The 
best scaling was found using K — 0.76,0.67,0.54,0.47 
for x = 0.1 . . . 0.4, respectively, and is shown in Fig.[4{a). 
The errors of K are discussed below. Despite the fact 
that K{x) was chosen to scale the T c domes, by the same 
token other critical temperatures scale as well. Fig. @Ja) 
also shows, T g /T™ ax and T N /T™ ax , for all families, as 
a function of Ap m . The T g /T™ ax curves of all families 
collapse onto each other and Tn jT™ ax for the x = 0.2 to 
0.4 families collapse onto each other. The reason Tjy of 
the x — 0.1 family is not in-line with the others is due to 
interactions between planes. 

The Neel temperature stems from the three dimen- 
sional interaction, and it is a function of not only the 
in-plane J, but also interplane couplings J± and other 
anisotropics. In the Heisenberg model Tn = Jt]y(a e ff) 
where the effective anisotropy ot e ff is mainly set by 
J_i_/J, and t jy ( a P j f ) is a known logarithmic function of 
the anisotropy [14[. We determine a e ff and extract J 
from T/v- This is done by measuring the temperature de- 
pendence of the muon rotation angular frequency ui(T), 
which is proportional to order parameter M(T). The 
function M(T/Tn)/Mq depends only on a e ff. The in- 
set of Fig. Eta) shows u(T/T n )/uj(T -> 0) for two sam- 
ples marked in Fig. Ufa), and a fit to the predicted be- 
havior given in Refs. [14j and The x — 0.1 sample 
clearly has a bigger a e ff than the x = 0.3 sample, and 
is more 3D like. Using this method we determined ct e ff 
for all samples with Neel order and defined the quan- 
tity Tjtf r = Tff /t^{a e ff) for these samples Q- For zero 
doping Tjf r = J . When the system is doped, Tjv is 
also affected by hopping and T^ r = J is no longer valid. 
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FIG. 4: (a) The CLBLCO phase diagram after rescaling: For 
each family (a;), the critical temperatures are normalized by 
T c at optimal doping, and y is replaced by Ap m (see text for 
details), (b) The same as (a) but the Neel temperature is 
corrected for anisotropy contribution so it is the same as J 
for the parent compound (see text), (c) The zero tempera- 
ture muon oscillation angular frequency as a function of Ap m 
for all four CLBLCO families; an equivalent to the staggered 
magnetization Mo versus mobile hole density. The arrows 
show the expected variation of the critical doping from a 5% 
variation in t/J. 



T^ r replaces T/v in Fig. HJb) which otherwise is the same 
as Fig. BJa). After this replacement the entire phase di- 
agram scales to a single unified curve, indicating that 
rpmax ^ j anc j a single energy scale controls mag- 
netism and superconductivity. 

Fig. |3Jc) shows u>(T — > 0) as a function of Ap m for 
each family. The scalability of the phase diagram, as 
explained above, suggests that Ap m is a parameter pro- 
portional to the mobile hole density variation. Hence 
Fig. 0Ja) is equivalent to a plot of the AFM order pa- 
rameter at zero temperature as a function of mobile hole 
density. This plot shows that the order parameter is uni- 
versal for all families, and in particular the AFM critical 



doping is family-independent. To demonstrate this point 
we show, using the two arrows in Fig. SJc), what should 
have been the difference in the critical doping had it been 
proportional to t/J, and changed between the x = 0.4 
and x — 0.1 by 5% (of 0.3) as indicated in Fig. OJc). 
Thus, we conclude that Mq(x, Ap m ) is x- independent, 
hence independent of t/J. 

The above conclusion could, a priori, depend on the 
choice of the if's and y m ax- A different set of K's or y m ax 
would shift the magnetic critical doping with respect to 
each other. However, it will also shift the normalized 
T c domes, T g line, and T^ r line with respect to each 
other. We have attempted to use a different set of K's 
and y m aa;'s, which will not noticeably destroy the scaling 
of the normalized critical temperatures. We could not 
generate a variation of more than 2% in the M critical 
doping. In other words, the different sets of K's and 
Umax's always kept the critical doping well between the 
two vertical arrows in Fig.[4jc). 

This surprising result could be discussed using three 
scenarios: (I) That different values of t/J correspond to 
different values of K, namely, the changes in t/J are can- 
celled out by the re-scaled doping axis. However, there is 
no theoretical backing for his scenario. (II) It may be that 
at low temperatures the effective Hamiltonian is given by 
a t-J model but with an effective t that is proportional 
to J. Indeed, there are indications that for small hole 
doping in an antifcrromagnet the bandwidth of the hole 
dispersion is set by J (instead of t) [lif . However, the 
spatial size of each hole quasiparticle (spin-polaron), and 
thus the critical doping, does depends on t/J. Moreover, 
it is not clear how spin-polarons destroy the AFM order 
and at which doping. (Ill) An alternative explanation is 
that the destruction of the AFM order parameter is not 
a result of single holes hopping and should be described 
by a completely different effective Hamiltonian; perhaps 
hopping of boson pairs fl6| . In this case t should be 
absorbed into the creation of tightly bound bosons leav- 
ing a prominent energy scale J. The proximity of the 
magnetic critical doping to superconductivity makes this 
possibility appealing. 

The last possibility could also solve a profound riddle 
in the study of the CLBLCO system. This system was 
found to obey the Uemura relation T c oc n s 17j, where 
n s is the superconducting carrier density, in both under- 
and overdoped regions 18| . At the same time T™ ax scales 
with J as indicated before. Therefore, the Uemura rela- 
tion should be rewritten as T c cx Jn s . What is then the 
role of i? Our finding that the magnetic order parame- 
ter versus doping is universal suggests that even before 
superconductivity appears, t becomes less relevant. This 
suggestion does have theoretical support . 

To conclude, an estimate of t and J from simple struc- 
ture considerations using neutron diffraction shows that 
the origin of the different energy scales between the 
CLBLCO families is mainly the different buckling angles. 



■5 



The difference in t/J between the two extreme families 
is about 5%. Although this is not an accurate way to 
measure the hopping rate or superexchange coupling, it 
does set the scale for the expected variation in the AFM 
critical doping. Using fiSR, the AFM order parameter as 
a function of oxygen was determined for different families 
of the CLBLCO system. We used a scaling transforma- 
tion to move from oxygenation level to mobile holes. Our 
measurements show that, at zero temperature, the order 
parameter as a function of mobile holes is independent 
of t/J within the required accuracy. 
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